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Abstract 

We consider the general physical situation of a quantum system T~Co interacting with a chain of 
exterior systems ®^V*7Y, one after the other, during a small interval of time H and following some 
Hamiltonian H on Hq ® T~l- We discuss the passage to the limit to continuous interactions (h — > 0) 
in a setup which allows to compute the limit of this Hamiltonian evolution in a single state space: a 
continuous field of exterior systems ® jp^+1~L. Surprisingly, the passage to the limit shows the necessity 
for 3 different time scales in H. The limit evolution equation is shown to spontaneously produce 
quantum noises terms: we obtain a quantum Langevin equation as limit of the Hamiltonian evolution. 
For the very first time, these quantum Langevin equations are obtained as the effective limit from 
repeated to continuous interactions and not only as a model. These results justify the usual quantum 
Langevin equations considered in continual quantum measurement or in quantum optics. We show that 
the three time scales correspond to the normal regime, the weak coupling limit and the low density 
limit. Our approach allows to consider these two physical limits altogether for the first time. Their 
combination produces an effective Hamiltonian on the small system, which had never been described 
before. We apply these results to give an Hamiltonian description of the von Neumann measurement. 
We also consider the approximation of continuous time quantum master equations by discrete time 
ones. In particular we show how any Lindblad generator is obtained as the limit of completely positive 
maps. 
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I. Introduction 

Quantum Langevin equations as a model for quantum open systems have been 
considered for at least 40 years (for example [FKM], [FLO], [AFL]). They have 
been given many different meanings in terms of several definitions of quantum 
noises or quantum Brownian motions (for example [G-Z], [H-P], [GSI]). One of the 
most developed and useful mathematical languages developed for that purpose is 
the quantum stochastic calculus of Hudson and Parthasarathy and their quan- 
tum stochastic differential equations ([H-P]). The quantum Langevin equations 
they allow to consider have been used very often to modelize typical situations of 
quantum open systems: continual quantum measurement ([Bal], [B-B]), quantum 
optics ([F-R], [FRS] [Ba2]), electronic transport [BRSW], thermalization ([M-R], 
[L-M]), etc. 

The justification for such quantum Langevin equation is often given in terms 
of some particular approximations of the true Hamiltonian interaction dynamic: 
rotating wave approximation, Markov approximation, large band approximation 
(cf [G-Z] chapter 11). 

They are also often justified as natural dilations of quantum master equations 
on the small system. That is, for any (good) semigroup of completely positive 
maps on the small system (with Lindblad generator £), one can dilate the small 
system with an appropriate Fock space, and obtain an explicit quantum stochastic 
differential equation on the whole space. The unique solution of this equation is a 
unitary evolution (in interaction picture) such that the trace on the small system 
of the induced evolution yields the original semigroup. This corresponds, at the 
quantum level, to the well-known way of realizing a concrete Markov process from 
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a given semigroup (or generator) by adding a noise space to the (classical) system 
space and solving an adequate stochastic differential equation. 

Some quantum stochastic differential equations have also been obtained in 
the so-called stochastic limit from explicit Hamiltonian dynamics ([A-L], [AGL], 
[ALV]). This shows some similarities with the results described here, but the limits 
considered in those articles are in the sense of the convergence of processes living 
in a different space than the one of the Hamiltonian dynamic. 

In this article we consider the effective Hamiltonian dynamic describing the 
repeated interactions, during short time intervals of length h, of a small system Tig 
with a chain of exterior systems ®iv*^- We embed all these chains as particular 
subspaces, attached to the parameter h, of a continuous field 

<g)W 

R+ 

in such a way that the subspaces associated to the chain increase and fill the 
field when h tends to 0. This framework may seem to specialize to the case of a 
zero-temperature exterior system; actually, as was noted by the first author and 
Maassen, it also applies to the case of positive temperature, using the cyclic (GNS) 
representation of the given state. 

By developing an appropriate language of the chain ®^v*7i and of the field 
^ffi+Ti, and by describing the discrete time Hamiltonan evolution generated by the 
repeated interactions, we are able to pass to the limit when h — > and prove that 
the limit evolution operator is the solution of a quantum stochastic differential 
equation. This limit is obtained in the strong topology of operators in a single 
space: the continuous field 7"^ an d implies the weak convergence of the 

Heisenberg evolutions of any observable. 

Of course, such a limit cannot be obtained without assumptions on the ele- 
mentary interaction Hamiltonian H. This is similar to the central limit theorem: 
a random walk gives a trivial limit when its time step h goes to zero and it is 
only when suitably renormalized (by a factor \fh) that it yields a Gaussian. Other 
normalizations give either trivial limits or no limit at all. 

In our Hamiltonian context the situation is going to be the same. For a non- 
trivial limit of these repeated interactions to exist, we will need the Hamiltonian H 
to satisfy some renormalization properties. The surprise here is that the necessary 
renormalization factor is not global, it is different following some parts of the 
Hamiltonian operator. We identify 3 different time scales in H: one of order 1, 
one of order y/h, one of order h. 

We describe a class of Hamiltonian which seems to be typical for the above 
conditions to be satisfied. These typical Hamiltonians are clearly a combination 
of free evolution, weak coupling limit typical hamiltonians and low density limit 
typical Hamiltonians. This physically explains the three different time scales. 
But the originality of our approach allows to consider both limits together; to 
our knowledge this constitutes a novelty in the literature. As a consequence, the 
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combination of the two limits shows an effective Hamiltonian for the small system 
which is very surprising: it contains a new term 

V* D~ 2 (sin D — D)V 

which comes from the presence of both the weak coupling and the low density 
limit in the Hamiltonian. It seems that such a term had never been described 
before. Again notice a possible extension of our results: only the case of time- 
independent coupling is discussed here but results for time-dependent ones can 
easily be deduced. 

This article is structured as follows. 

In section II we present the exact mathematical model of repeated quantum 
interactions and end up with the associated evolution equation (subsection II. 1). 
We then introduce a mathematical setup for the study of the space (S)iv* ^ which 
will help much for passing to the continuous field. In particular this includes a 
particular choice of an orthonormal basis of the phase space and a particular choice 
of a basis for the operators on that phase space (subsection II. 2). Finally we show 
how the typical evolution equations obtained in II. 1 are the general model for 
the unitary dilation of any given discrete semigroup of completely positive maps 
(subsection II. 3). 

Section III is devoted to presenting the whole formalism of the continuous 
atom field. In subsection III.l we present the space which is candidate for rep- 
resenting the continuous field limit of the atom chain. It is actually a particular 
Fock space on which we develop an unusual structure which clearly shows the 
required properties. In subsection III. 2 we present the natural quantum noises on 
the continuous field and the associated quantum stochastic integrals, the quantum 
Ito formula and the quantum stochastic differential equations. In subsection III. 3 
we concretely realize the atom chain of section II as a strict subspace of the atom 
field. Not only do we realize it as a subspace, but also realize the action of its 
basic operators inside the atom field. All these atom chain subspaces are related 
to a partition of M + . When the diameter of the partition goes to 0, we show 
that the corresponding subspace completely fills the continuous field and the basic 
operators of the chain converge to the quantum noises of the field (with convenient 
normalizations). Finally, considering the projection of the continuous atom field 
onto an atom chain subspace, we state a formula for the projection of a general 
quantum stochastic integral. 

In section IV all the pieces of the puzzle fit together. By computing the 
projection on the atom chain of a quantum stochastic differential equation we 
show that the typical evolution equation of repeated interactions converges in 
the field space to the solution of a quantum Langevin equation, assuming the fact 
that the associated Hamiltonian satisfies some particular renormalization property 
corresponding to three different time scales. It is to that result and to some of its 
extensions that subsection IV. 1 is devoted. In subsection IV. 2 we describe a family 
of Hamiltonians which seems to be typical of the conditions obtained above. We 
show that this family of Hamiltonians describes altogether free evolution, weak 
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coupling limit and low density limit terms. Computing the associated quantum 
Langevin equation at the limit, we obtain an effective Hamiltonian on Hq which 
contains a new term. This new term appears only when weak coupling and low 
density limits are in presence together in the Hamiltonian. In subsection IV. 3, 
we apply these results to describe the von Neumann measurement apparatus in 
the Hamiltonian framework of repeated quantum interactions. In subsection IV. 4 
we explicitly compute a simple example. In subsection IV. 5 we show that our 
approximation theorem puts into evidence a natural way that completely positive 
maps have to converge to Lindblad generators. 

II. Discrete dynamics on the atom chain 

II. 1 Repeated quantum interactions 

We here give a precise description of our physical model: repeated quantum 
interactions. 

We consider a small quantum system Hq and another quantum system H 
which represents a piece of environment, a measuring apparatus or incoming 
photons. . . We consider the space Hq <8> H in order to couple the two systems, 
an Hamiltonian FT on Hq <g> H which describes the interaction and the associated 
unitary evolution during the interval [0, h] of time: 

TL = e~ ihH . 

This single interaction is therefore described in the Schrodinger picture by 

p I— > TL pTL* 

and in the Heisenberg picture by 

X ^ TL*XTL. 

Now, after this first interaction, we repeat it but this time coupling the same 
Hq with a new copy of H. This means that that new copy was kept isolated 
until then; similarly the previously considered copy of H will remain isolated for 
the rest of the experience. One can think of many physical examples where this 
situations arises: in repeated quantum measurement where a family of identical 
measurement devices are presented one after the other before the system (or a 
single device is refreshed after every use), in quantum optics where a sequence 
of independent atoms arrives one after the other to interact with a field in some 
cavity for a short time. More generally it can be seen as a good model if it is 
assumed that perturbations in H due to the interaction are dissipated after every 
time h. 

The sequence of interactions can be described in the following way: the state 
space for the whole system is 

Hq®®H 

]N* 
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Index for a few lines only the copies of H as Hi, H2, ■ ■ • Define then a unitary 
operator Ej n as the canonical ampliation to Ho <S> Hi <S> H2 ® ... of the operator 
which acts as iL on Ho ®H n ; that is, TL n acts as the identity on copies of H other 
than H n . 

The effect of the n-th interaction in the Schrodinger picture writes then 

p ^ Ij n plj* n , 

for every density matrix p, so that the effect of the n first interactions is 

p h-> u n pu* n 

where (u n ) ne ]N is a sequence in B(Ho <E> (SW* ^0 which satisfies the equations 

J u n+l = ILn+l U n s 

I uo = I. [L> 
It is evolution equations such as (1) that we are going to study in this article. 

II. 2 Structure of the atom chain 

We here describe some useful mathematical structure on the space ®jn*H 
which will constitute the main ingredient of our approach. 

Let us fix a particular Hilbertian basis (X^gAufo} f° r the Hilbert space H, 
where we assume (for notational purposes) that ^ A. This particular choice of 
notations is motivated by physical interpretations: indeed, we see the X 1 , i e A, 
as representing for example the different possible excited states of an atom. The 
vector X° represents the "ground state" or "vacuum state" of the atom and will 
usually be denoted O. 

Let T$ be the tensor product (S>iv* ^ with respect to the stabilizing sequence 
O. In other words, this means simply that an orthonormal basis of T£> is given by 
the family 

{X A - AeVjn^a} 

where 

- the set V]n,a is the set of finite subsets 

{(m,ii),...,(rifc,ifc)} 

of IN* x A such that the n^s are mutually different. Another way to describe the 
set ?jv*,a is to identify it to the set of sequences {A n ) ne]N , with values in AU {0} 
which take a value different from only finitely often. 

- Xa denotes the vector 

O <g> . . . <g> O <g> X h <g) Q <g) . . . <g) Q <g) X i2 (g) ... 

where X %1 appears in ni-th copy of H... 

The physical signification of this basis is easy to understand: we have a chain 
of atoms, indexed by IN*. The space T$ is the state space of this chain, the 
vector Xa with A = {(ni, ii), . . . , (n^, ik)} representing the state in which exactly 
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k atoms are excited: atom m in the state %\, etc, all other atoms being in the 
ground state. 

This particular choice of a basis gives T$ a particular structure. If we denote 
by Tl> n ] the space generated by the Xa such that Ac. {1, . . . , n} x A and by T$[ m 
the one generated by the Xa such that A C {m, m + 1, . . .} x A, we get an obvious 
natural isomorphism between T$ and T$ n _i] ® T3>[ n given by 

[/^P) = /(An{l,...,n-l}xA) ff(An{n,...}xA). 

Put {a*; i,j e A U {0}} to be the natural basis of 13(H), that is, 

a}(X*) = <y ifc X''. 

We denote by o*(n) the natural ampliation of the operator a* to T$ which acts 
on the copy number n as a* and the identity elsewhere. That is, in terms of the 
basis Xa, 

aj(n)X A = \n,i)eAX A\{n,i)U(n,j) 

if neither i nor j is zero, and 

al(n)XA = \n,i)£AXA\(n,i), 
a®(n)XA = l(n,0)6Au(nj), 
a o( n )XA = \n,Q)&A^Ai 

where (n, 0) G A actually means "for any i in A, (n, i) A" . 

II. 3 Unitary dilation of completely positive maps 

The evolution equations 

u n = Ej n . . . iLi 

obtained in the physical setup of repeated quantum interactions are actually of 
mathematical interest on their own for they provide a canonical way of dilating 
discrete semigroups of completely positive maps into unitary automorphisms. 

The mathematical setup is the same. Let iL be any operator on Ho <g> H. Let 
T$> = ®n*H and (Ej n ) ne]N * be defined as in the above section. We then consider 
the associated evolution equations 

u n = lL n ...lL 1 (1) 

with uq = I. 

The following result is obvious. 

Proposition 1.- The solution (a n ) ne]N of (1) is made of unitary (resp. isometric, 
contractive) operators if and only if IL is unitary (resp. isometric, contractive). 

■ 

Note that if iL is unitary, then the mappings 

j n {H) = u* n Hu n 



7 



Stephane ATTAL and Yan PAUTRAT 



are automorphisms of B(Ho <g) H). 

Let 1E be the partial trace on Ho defined by 

for all 4>,ijj E TCo and every operator i7 on Ho <S> T3>. 

Unitary dilations of completely positive semigroups are obtained in the fol- 
lowing theorem. Recall that, by Kraus' theorem, any completely positive operator 
£ on B(Hq) is of the form 

£(X) = A*XA> 

where the summation ranges over (Au{0}) 2 , the Ai are bounded operators and the 
sum is strongly convergent. Conversely, any such operator is completely positive. 
Remark: Of course the Kraus form of an operator is a priori indifferent to the 
specificity of the value i = 0. The special role played by one of the indices will 
appear later on. 

Theorem 2. - Let E, be any unitary operator on Ho®H. Consider the coefficients 
(iL*)j j 6 Au{o}; which are operators on TCo, of the matrix representation of Ej in the 
basis Q,X\ ieAofH. 

Then, for any X e 13 (Ho) we have 

E [j n (X®I)]=£ n (X) 

where £ is the completely positive map on B(Hq) given by 

£(X)= 

i£AU{0} 

Conversely, consider any completely positive map 

£(X)= A *i XA i 

i£AU{0} 

on B(Ho) such that £(I) = I. Then there exists a unitary operator Ej on Ho <8> H 
such that the associated unitary family of automorphisms 

jn(H) = u* n Hu n 

satisfies 

E [jn{X®I)]=£ n {X), 

for all neW. 
Proof 

Consider E = (iL*))j jeAU { } such as in the above statements. Consider the 
unitary family 

u n = E n . . .Ei. 

Note that 

Un+l = E n+ iU n . 
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2L° TL\ ... 

V ! i '••/ 



Put j n (H) = u* n Hu n for every operator H on Ho <8> H- Then, for any operator X 
on Ho we have 

j n +i(X <g> I) = u* n E* n+1 (X <g> I)E n+1 u n . 

When considered as a matrix of operators on Ho, in the basis fi, X\ i e A of H, 
the matrix associated to X <g> I is of diagonal form. We get 

® /)&n+l = 

K)* (^!)* .-A /x o 

(ij)* (^)* ... X 

which is the matrix L n+1 (X) = (BUX)). . with 

5j(X) = £ (HTxili 

keAU{0} 

Note that the operator L n+1 (X) acts non trivialy only on the tensor product of 
Ho with the (n + l)-th copy of H. When represented as an operator on 

H <g> T$ n+1] = (Ho ® T£ n] ) <g> W 

as a matrix with coefficients in £>(7Y <8> T3> n ]) it writes exactly in the same way as 
above, just replacing (which belongs to B(Ho)) by 

B«(X)®J| TOa] . 

Also note that, as can be proved by an easy induction, the operator u n acts 
on Ho <8> T$ n ] only. As an operator on Ho <8> T$ n+1 ] it is represented by a diagonal 
matrix. Thus j n+ i(X) = u^Li n +i(X)u n can be written on Ho <8> T£>( n+ i)] = 
7i <8> T$ n ] ® H as a matrix of operators on Ho <8> T$ n ] by 

n+1 (X®I));=j n (Bj(X)®7). 

Note that 5g(X) = ^ . 6Au{0} (iL - )*XiL° which is the mapping £(X) of the state- 
ment. 

Put T n (X) = IE [j n (X <g> /)]. We have, for all 0, * e 7Y 

= < ® o , (jn(s;.(A:) ® /))..* ® ft > 

= <0®nii. n] ®o w , (j n (sj(x)® i)).^*® ®o w > 

= < $ ® f2 TOn] , j n (^(x) ® /)* ® n TOn] > 

= <0,T n (£pT))#>. 

This proves that T n+1 (X) = T n (£(X)) and the first part of the theorem is proved. 
Conversely, consider a decomposition of a completely positive map £ of the 

form 



C(X)= Yl A * XA i 



i£AU{0} 
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for a familly (^4f)ieAu{o} °f bounded operators on 7io such that 

i£AU{0} 

We claim that there exists a unitary operator iL on 7Yo <8> W of the form 

A> 

iL = I 



Indeed, the condition X]ieAu{o} A\Ai = I guarantees that the first columns of iL 
are made of orthonormal vectors of T~Lq ® 7Y. We can thus complete the matrix by 
completing it into an orthonormal basis of Ho^Ti,- This makes out a unitary matrix 
iL the coefficients of which we denote by (Aj). j €Au ^ y Note that A® = A i+1 . We 
now conclude easily by the first part of the theorem. 



Ill From the atom chain to the atom field 

III.l Structure of the atom field 

We now describe the structure of the continuous version of Ti>. The struc- 
ture we are going to present here is rather original and not much expanded in the 
literature. It is very different from the usual presentation of quantum stochastic 
calculus ([H-P]), but it actually constitutes a very natural language for our pur- 
pose: approximation of the atom field by atom chains. This approach is taken 
from [Atl]. 

We first start with a heuristical discussion. 

By a continuous version of the atom chain T$ we mean a Hilbert space with 
a structure which makes it the space 



<S> = Q<)H. 

R+ 

We have to give a meaning to the above notation. This could be achieved by in- 
voquing the framework of continous tensor products of Hilbert spaces (see [Gui]), 
but we prefer to give a self-contained presentation which fits better with our ap- 
proximation procedure. 

Let us make out an idea of what it should look like by mimicking, in a con- 
tinuous time version, what we have described in T£>. 

The countable orthonormal basis Xa, A e Vjn*,a is replaced by a continuous 
orthonormal basis dxa, & £ Vjr,K: where Vjr^a is the set of finite subsets of M + x A. 
With the same idea as for T$, this means that each copy of 7i is equipped with an 
orthonormal basis O, d\\, « G A (where t is the parameter attached to the copy we 
are looking at). The orthonormal basis above is just the one obtained by specifying 
a finite number of sites ti, . . . , t n which are going to be excited, the other ones 
being supposed to be in the fundamental state O, and by specifying their level of 
excitation. 
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The representation of an element / of T£>: 

\\f\\ 2 = E i/(^)i 2 

is replaced by an integral version of it in $: 

\\f\\ 2 =[ \ffda. 

This last integral has to be explained: the measure da is a "Lebesgue measure" 
on V]R t A, as will be explained later. From now on, the notation V will denote, 
depending on the context, spaces of the type Pjv*,a or Vu,k- 

A good basis of operators acting on $ can be obtained by mimicking the 
operators a*(n) of T$. We will here have a set of infinitesimal operators da l j(t), 
i,j e A U {0}, acting on the "t-th" copy of H by: 

da®(t) dxa = dxa dt t t ga 

da°(t) dxa = dXau{(t,i)} %£o 

da l (t) dxa = dxa\{(t,i)} dt \t,i)z<j 

da)(t) dXa = dXa\{(t,i)}U{(t,j)} \t,i)£a 

for all i,j E A. 

We shall now describe a rigourous setup for the above heuristic discussion. 

We recall the structure of the bosonic Fock space $ and its basic structure 
(cf [Atl] for more details and [At3] for a complete study of the theory and its 
connections with classical stochastic processes). 

Let TL be, as before, a Hilbert space with an orthonormal basis A 1 , % e Au{0} 
and let H' be the closed subspace generated by vectors A 1 , j G A (or simply said, 
the orthogonal of A ). 

Let $ = r s (L 2 (iR + , Tl')) be the symmetric (bosonic) Fock space over the 
space L 2 (IR ,7i'). We shall here give a very efficient presentation of that space, 
the so-called Guichardet interpretation of the Fock space. 

Let V (= Vm,a) be the set of finite subsets {(si, h), ■ ■ . , (s n) in)} of M + x A 
such that the Sj are two by two different. Then V = U„?„ where V n is the 
subset of V made of n-elements subsets of 1R + x A. By ordering the iR + -part 
of the elements of a G V n , the set V n can be identified to the increasing simplex 
E n = {0 < ti < ■ ■ ■ < t n } x A of M n x A. Thus V n inherits a measured space 
structure from the Lebesgue measure on M n times the counting measure on A. 
This also gives a measure structure on V if we specify that on Vo = {0} we put the 
measure 8$. Elements of V are often denoted by a, the measure on V is denoted 
do. The a-field obtained this way on V is denoted T . 
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We identify any element a G V with a family {<Tj, % G A} of (two by two 
disjoint) subsets of M + where 

Oi = {s G iR + ; (s,i) G a}. 

The Fock space $ is the space L 2 (V, da). An element / of <3> is thus a 
measurable function / : V — > C such that 

||/|| 2 = / \f(a)\ 2 da <oc. 

One can define, in the same way, V[ a ,b] and <£•[«,&] by replacing iR + with [a, 6] C iR + . 
As in discrete time, there is a natural isomorphism between $[o,t] <8> $[t,+oo[ given 
hy h® f where /(a) = n [0, t])#(cr n (t, +oo[). 
We shall use the following notations: 

®t]=®[0,t], $[t = ®[t,+oo[- 

Define O to be the vacuum vector, that is, Q(a) = 60(a). 

We now define a particular family of curves in $, which is going to be of great 
importance here. Define x£ G $ by 

x \a) = ( %*](*) if a = 

1 otherwise. 

Then notice that for all t G M + we have that xl belongs to $[o,t]- We actually 
have much more than that: we have 

Xl - Xl e $[s,t] for all s < t. 
This last property can be checked immediately from the definitions, and it is 
going to be of great importance in our construction. Also notice that xl an d X 3 S 
are orthogonal elements of $ as soon as i 7^ j. As we will see later on, apart from 
trivialities, the curves (Xt)t>o are ^ ne om y ones to share these properties. 

These properties allow to define the so-called Ito integral on $. Indeed, let 
g = {(gl)t>o, i G A} be families of elements of $ indexed by both M + and A, 
such that 

i) 1 1— > ||<7||| is measurable, for all i, 

ii) 0je$[ O ,t] for all t, 

then one says that (7 is Ito integrable and we define its Ito integral 



ieA 

to be the limit in $ of 



J' 



where 5 = {tj, j'G-SV} is a partition of iR + which is understood to be refining 



and to have its diameter tending to 0, and (g l .)i is an Ito integrable family in $, 
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such that for each z, t t— > gl is a step process, and which converges to (g^i in 
L 2 (M + xV). 

Note that by assumption we always have that g l t . belongs to <& tj .j and Xt j+1 ~ 
Xtj belongs to Qfajj+i], hence the tensor product symbol in (2). 
Also note that, as an example, one can take 



mi?) 



gl = -, r / P tj g l 8 ds 

where Pt is the orthogonal projection onto $[o,t]- 

One then obtains the following properties ([Atl], Proposition 1.4): 

Theorem 3.- The Ito integral 1(g) = £\ J °° gl d\\, of an Ito integrable family 
g = {g l )i£\, is the element of Q given by 

if a = 

£va( CT \ (V<M)) ifVo-ea l . 
It satisfies the Ito isometry formula; 

iia<7)ii 2 = E/ ^ = £/ Ml dt - ( 3 ) 



In particular, consider a family / = (f l ) ie \ which belongs to L 2 (Vi) = 
L 2 (M + x A), then the family (/'(t)O), t G M + , i G A, is clearly Ito integrable. 
Computing its Ito integral we find that 



poo 

J (/) = E/ d X j 



i6A ' 

is the element of the first particle space of the Fock space $ associated to the 
function /, that is, 

ifa = {s}, 
otherwise. 



'(/)(*) 



Let us define the "adjoint" mapping of the Ito integral. For all / G $, all i in 
A and all t G M + , consider the following mapping on P: 

[D*f] (o-) = f(*U{( S ,i)})l aC[0A . 

We then have the following result ([Atl], Theorem 1.6). 

Theorem 4. [Fock space predictable representation property] -For all f G 
all i E A and for almost all t G IR + , the mapping D l t f belongs to $ = L 2 (V). 
Furthermore, the family (D l ,f)i is always Ito integrable and we have the represen- 
tation 

/ = /(0)n + E/ D if*x\ ( 4 ) 

zeA- 70 
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with the isometry formula 

poo 

" = \fm\ 2 + J2 \\ D if\\ dt - (5) 



As an immediate corollary we get the following. 

Corollary 5.- The representation (4) of f is unique. In particular, if g G $ is of 
the form 



ieA 

then for almost all t, all i in A, 



/>oo 

g = cn + V / h\d x l 

iEA J ° 

Dig = K • 

Let f<EL 2 (V n ), one can easily define the iterated Ito integral on $: 



In(f)= [ f(0-)d X a 

JV-r, 



'v. 

by iterating the definition of the Ito integral: 



J n(f)= E / / •••/ f 1 '-' in (tu---,tn)ndxi\ ...dxt 

n,...,i n €A J o Jo Jo 

We obtain this way an element of $ which is actually the representant of / in the 
n-particle subspace of $, that is 

[I n (f)]( a ) = I /" '" (*i> •••,*») ^ o- = U . . . U {t n },„ 

\ otherwise. 

For any f £ V we put 

/ /(*)dXa 

to denote the series of iterated Ito integrals 

»=li ll ...,i„€A" / ° J ° J ° 

We then have the following representation ([Atl], Theorem 1.7). 

Theorem 6. [Fock space chaotic representation property] - Any element f of $ 
admits an abstract chaotic representation 



/ = / f(<r) d X a (6) 

JV 

satisfying the isometry formula 

r \f(o-)\ 2 da (7). 



V 
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This representation is unique. 

■ 

The above theorem is the exact expression of the heuristics we wanted in order 
to describe the space 

$ = (g)W. 

R+ 

Indeed, we have, for each t G iR + , a family of elementary orthonormal elements fi, 
dx\, i G A (a basis of Ti) whose (tensor) products d\ a form a continuous basis of 
$ (formula (6)) and, even more, form an orthonormal continuous basis (formula 
(7))- 

The attentive reader will have noticed that the only property of the curves 
X* that we really used is the fact that xl ~ Xl belongs to $[ s ,t] for all s < t. 
One can naturally wonder if there exists another such family, which will then 
allow another Ito integral and furnish another continuous basis for $ via another 
chaotic expansion property. 

Of course there are obvious curves that can be obtained from the x % \ for any 
function / on 1R + and any g G L 2 (1R + x A) put 

for all t G M + . Then one easily checks that (y.) satisfies the same property, 
namely, yt — y s belongs to &[ s ,t] f° r au s <t. But clearly the Ito integration theory 
obtained from y is the same as the one from x, except that scalar factors g l (s) will 
appear in the integration. 

One can wonder if there exist more complicated examples, giving rise to a 
different Ito integration. The following result shows that there are no more exam- 
ples. In particular, there is only one Ito integral, one chaotic expansion and one 
natural continuous basis ([Atl], Theorem 1.8). 

Theorem 7.- Let (yt) t > be a curve in $ such that y t — y s belongs to &[ s ,t] f or 
all s < t. Then there exist a function f on M + and g G L 2 (1R + x A) such that 

yt = f(t)n + Yl /V(s)fW s 
i J ° 

for allt G M + . 

■ 

III. 2 The quantum noises 

The space $ we have constructed is the natural space for defining quantum 
noises. These quantum noises are the natural, continuous-time, extensions of the 
basis operators a* (n) we met in the atom chain T$. 

As indicated in the heuristic discussion above, we shall deal with a family of 
infinitesimal operators da l Jt) on $ which act on the continuous basis dx a m the 



15 



Stephane ATTAL and Yan PAUTRAT 



same way as they discrete-time counterparts a l Jn) act on the Xa- The integrated 
version of the above heuristic infinitesimal formulas easily gives an exact formula 
for the action of the operators a l -{t) on $: 

[«?(*)/](*)= £/(*\W«), 

[«o(*)/K*) = tfiWith) *»> 
Jo 

[<#)/>)=£ /(^\WiUW0 

s<i 

[a°(t)/](<7) = t/» 

for i, j ^ 0. 

All these operators, except a(j(t), are unbounded. But note that a good 
common domain to all these operators is 

J \<t\ |/(a)| 2 da<oo} . 

This family of operators is characteristic and universal in a sense which is close 
to the one of the curves \\. Indeed, one can easily check that in the decomposition 
of $ ~ $ s j $[s,t] ® &[t, the operators a*(£) — a*-(s) are all of the form 

7(8)(aj(t)-a5(a))|* [M] (8)/. 

This property is fundamental for the definition of the quantum stochastic integrals 
and, in the same way as for (%?), these operator families are the only ones to share 
that property (cf [Coq]). 

This property allows to consider Riemann sums: 

£ff tfc (aj(f fc+ i)-aj(f*)) (8) 
k 

where S = {0 = to < t\ < . . . < t k < . . .} is a partition of M + , where (H t ) t>Q is a 
family of operators on $ such that 

- each H t is an operator of the form H t <8> I in the tensor product space 
$ = (g) (we say that Jit is a t-adapted operator and that (H t ) t>0 is an 
adapted process of operators), 

- (H t ) t>0 is a step process, that is, it is constant on intervals: 

k 

and where the operator product H tk (aUtk+i) — a l j(tk)) is actually a tensor prod- 
uct of operators 

H tk ®{a)(t k+1 )-a)(t k )). 
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Note that, in particular, the above "product" is commutative and does not impose 
any new domain constraint. 

The resulting operator associated to the Riemann sum (8) is denoted by 

f oo 

Hsda'is). 

If we denote by T the above operator and by T t the operator 

/•OO 

H s da)(s) = / H s t m (s)da)(s) 
Jo 

we can then compute the action of T on a "good" vector / of its domain and we 
obtain (cf [A-M] for more details) 

poo poo 

Tf = J2 T t D?fd X 1+ H t Dlfdxl (9) 
keA Jo Jo 

with the notations: D° = P t and dxt = dt. For general operator processes (still 
adapted but not step process anymore) and general /, it is equation (9) which is 
kept as a definition for the domain and for the action of the operator 



/•OO 

T= / HsdaUs). 
Jo 



The maximal domain and the explicit action of the above operator can be described 
but is not worth developing here. The interested reader may refer to [ At3] , chapter 
12 or to [A-L]. There are particular domains where the definition simplifies. The 
one we shall use here is the case of coherent vectors. 

Indeed, if 4> is an Y element of L 2 (M + ,H'), consider the associated coherent 
vector e((p) in That is, 

=n 

Put <j>o(s) = 1 for all s. If (p is such that 

/ |</> J (s)| (1 ~ <5oi)(2 ~ <5 ° 3) || #,£(</)) || 2 " 5oj ds <C oo 
Jo 

then Jq H s da l j(s) is well-defined on e{4>) with 

<e(ip), [ H s da)(s)e((j))>= [ ^ J (s)(p l (s)< e{ip) , H s e{4>) > ds. 
Jo ' Jo 

for all i) G L 2 (R + ,W). 

III. 3 Embedding and approximation by the Toy Fock space 

We now describe the way the atom chain and its basic operators can be 
realized as a subspace of the Fock space and a projection of the quantum noises. 
The subspace associated to the atom chain is attached to the choice of some 
partition of M + in such a way that the expected properties are satisfied: 

- the associated subspaces increase when the partition refines and they con- 
stitute an approximation of $ when the diameter of the partition goes to 0, 
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- the associated basic operators are restrictions of the others when the parti- 
tion increases and they constitute an approximation of the quantum noises when 
the diameter of the partition goes to 0. 

Note that this approximation has deep interpretations in terms of approxi- 
mations of n-dimensional classical noises by extremal random walks in M n whose 
jumps take n + 1 different values. This aspect is developed in [A-P]. 

Let S = {0 = to < t\ < • ■ ■ < t n < ■ ■ ■} be a partition of M + and 
5(S) = supjti+i — ti\ be the diameter of S. For S being fixed, define <E> n = 
^[t„_i,t„]> ndlN* . We clearly have that $ is naturally isomorphic to the countable 
tensor product ® n6 jv*^ ) n (which is understood to be defined with respect to the 
stabilizing sequence (0) n6 ^v). 

For all nGlV*, define for i,jeA 

xl — yl 
X\n) = " _ tn ~ 1 e $ n , 

y^n tn — 1 

4(n) = ^M^tiloP 1] , 
a){n) = Pi] o (a)(t n ) - a}(* n _i)) o P 1} , 
a°(n) = P 1] o^^fei) , 

y*n tn— 1 

ag(ra) = P ], 

where for i = 0, 1, Pj] is the orthogonal projection onto L 2 {Vi) and where the above 
definitions are understood to be valid on $ n only, the corresponding operator 
acting as the identity operator / on the others $ m 's. 

For every A e V = Pjv*,a, define from the X l (n)'s in the same way as 
for T$: 

x A = tt <g> . . . <g> O <g> JT 1 (m) <g> o <g> . . . <g> O <g> X i2 (n 2 ) <g> . . . 

in ®n6JV*^n- 

Define T$(iS) to be the space of /e<& which are of the form 

Aer 

(note that the condition ||/|| 2 = YIaev \ f(A)\ 2 < oo is automatically satisfied). 

The space T$(iS) is thus clearly identifiable to the spin chain H>. The space 
T<E>(<S) is a closed subspace of <E>. We denote by Ps the operator of orthogonal 
projection from $ onto T<3>(<S). One can prove for example that the projection of 
an exponential vector is an "exponential vector" of the embedded toy Fock space: 
indeed, a direct computation shows that for any 4> in L 2 (M + ,TC'), 

(Pse(J>))(A) = l[ J] Mn) 

i nEAi 
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where the function <fi belongs to I 2 (IN*, H') and is defined by 

1 f tn 
<t>i( n ) = —jf . I 4>i(s)ds. 

ytn tn-1 Jt n -i 

We will denote by e(4>) such a discrete time version of a coherent vector. 

We shall now check that the above operators a* (n) act on r H>( l S) in the same 
way as the the basic operators of T£>. 

Proposition 8. - We have, for all i,jeA 

a\{n) Xi{n) = <%0 

a i n = o 

a)(n)X k (n)=5, k X3(n) 
a} O = 



a°(n) Xi{n) = 
af(n)n = X^n) 

a° 9 (n)X k (n) = 
a) = 0. 

Proof 

This is a direct application of the definitions and computations using equation 
(9), cf [At 2] for details. For example: 

ai(n)Xi(n) = —— (a£(t n ) - a^U-x)) / O d X { 

1 



1 



rtn ftn 

/ H(t) - 4(tn-i)) n d x k t + / SijQ dt 

Ur- A J tn — 1 Jin— 1 



(0 + (t n - t n -i)SijQ) = SijQ. 



t n t 

And so on for the other cases. ■ 

Thus the action of the operators a* on the X % (n) is exactly the same as the 
action of the corresponding operators on the spin chain of section II; the operators 
a l j(n) act on T$(5) exactly in the same way as the corresponding operators do on 
T£>. We have completely embedded the toy Fock space in the Fock space. 

The action of operators a* (n) on discrete exponential vectors as defined above 
will be most useful in the sequel. The following lemma is deduced immediately 
from Proposition 8. 

Lemma 9. -For any </>, tp in L 2 (iR + , H) and for any t n -adapted operator H n the 
bracket 

'e(<j>), # n a;.(n+l)e(V>y 
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is equal to 

4>j(n + l)4>i(n + 1) (e(<j> n] ), H n e(^ n] )^ (e(4> [n+2 ), e(^ [n+2 )^ . 

This lemma is the basis for our future computations involving discrete-time quan- 
tum stochastic integrals (for more precise treatment of this subject see [Pal] or 
[Pa3]). 

We are now going to see that the Fock space $ and its basic operators a^(£), 
i,j G A U {0} can be approached by the toy Fock spaces T$(<S) and their basic 
operators a* (n). 

We are given a sequence (<S n ) n€ jv of partitions which are getting finer and finer 
and whose diameter 5(S n ) tends to when n tends to +oo. Let T$(n) = T<E>(<S n ) 
and P n = Ps n , for all neW. 

Theorem 10. 

i) The orthogonal projectors P n strongly converge to the identity operator I 
on <E>. That is, any f G $ can be approached in $ by a sequence (f n ) n€ w su °h 
that f n G T$(n) for aline IN. 

ii) If S n = {0 = t% < t\ < ■ ■ ■ < q < ■ ■ ■}, then for all teM + , the op- 
erators E fc . t n< t a}(A0, Y, k; t»<t y/tk-t^aiik), E k -t-<t V f k ~ f k-i a Uk), and 
^2k-t n <t( t7 k ~~ tfi-i) a o(^) converge strongly on V to a l j(t), a l (t), a°(£) and a°(f) 
respectively. 

Proof 

i) As the S n are refining then the (P n ) n forms an increasing family of orthog- 
onal projection in $. Let P^ = V n P n . Clearly, for all s < t, we have that xl ~ Xl 
belongs to RanPoo. But by the construction of the Ito integral and by Theorem 
6, we have that the xl ~~ Xl generate $. Thus P^ = I. 

ii) Let us check the case of a°. A direct computation shows that, for feV 
E yJtk-tk-i°i(Qf\(<r)= E 1 k«n[t S _ 1 ,t S ]|=i E /(^\Wi)- 

Put t n = sup{t£G<S n ;q < t}. We have 



E yJW -afrt) \f\\ 

k;tl<t J 

= | E 1 kin[*j_ 1 ,tj]i=i E f( a \i s }i)- E f( a \i s }i) 

<2 I I f(<r\{s}i)\ 2 d* + 2 j | E 1 k,n[t S _ 1> * J ]|>2 

x E /(*\wo 
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For any fixed a, the terms inside each of the integrals above converge to when n 
tends to +00. Furthermore we have, for n large enough, 



secrn[t n ,t] 



[ : E /(*\W0 2 *r< / W\ E l/(*\W0| 2 d* 

Jv 1 ,1 JV serT 

s<i + l 

= [ t+1 [ (\a\ + l)\f(a)\ 2 dads 
Jo Jv 

<(t+l) f (\a\ + l)\f(a)\ 2 da 
Jv 



which is finite for /GX>; 

E 1 kin[t«_ 1 ,t«]|>2 
k-ti<t 



E /(*\Wi)| da 



< 



< 



j ( E a k,n[^_ 1 ,^]i>2 1 E /( f7 \wo|) 

V k-,t»<t se<r i n[t™_ 1 ,t%] 

E E i/(*\woi) 2 d* 

k;t?<t s€<r;n[t£_ 15 t£] 

E i/(^\w*)i) 2 ^ 



da 



s<t n 



[ w\ E i/(^\Wi)i 2 ^ 

Jv ,^„. 



3<t n 



<(t+l) f (|a| + l)|/(a)| 2 rfa<cx) 
Jv 



in the same way as above. So we can apply Legesgue's theorem. This proves the 
result. 

The other cases are treated in the same way. See [At2] for details. ■ 

We have fulfilled our duties: not only the space T$(<S) recreates T£> and its 
basic operators as a subspace of $ and a projection of its quantum noises, but, 
when 5(S) tends to 0, this realisation constitutes an approximation of the space 
$ and of its quantum noises. 

To any operator H on $ we can associate the projected operator P$HPs 
which acts on the atom chain only and which approximates H (if H is bounded 
for example). 

We wish to compute the corresponding projections of the quantum stochastic 
integral operators. We reduce our computations to the case where integrals are 
of the type H = J °° R\ daUs), with 7^ (0,0), and satisfy the following 

conditions (HS): 

- the operator H is bounded and 
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- the integrands H]{t) are bounded for all t and t ^ \\Hj(t)\\ is square 
integrable if one of % or j is zero, essentially bounded otherwise. 

Even though they are rather restrictive, these hypotheses will suffice for our 
neeeds. 

The following result is a consequence of the theory and the computations 
developed in [Pal] and [Pa2] (adapted here to the case of higher multiplicity) and 
we do not reproduce the proof here. It is also stated in the following form in [Pa3] , 
chapter 4. 

Theorem 11. Let (z, j) 7^ (0,0) be fixed. Let H = J Hj(t)daUt) be a quantum 
stochastic integral on $ that satisfies the assumptions (HS). Then PsHPs is an 
operator on T$> of the form 

£5>f(n)a?(n+l) 

k,l n 

where the sum is over all couples (k,l) in (A U {0}) 2 different from (0,0) and is 
meaningful in the weak sense. The operators hf are given by : 

- if both i and j are nonzero, 

1 rt n +i 
hf(n) = 8 ki 5 lr -P s / Pt n H)(t) dt 

-i/i = 0, 

1 ftn+1 

hi (n) = Sij ,. . Ps / Pt n H°(t) dt 

V l n+1 — In Jt n 

and for all k 7^ 0, 

1 ftn+l 

hf(n) = 5 lr -P S / P tn H°(t) (a° k (t) - a° k (t n )) dt 



if 3 = 0, 

h k (n) = 5 ki - 1 P s i Pt n H l {t) dt 



tn+l 



and for all I 7^ 0, 



1 rt n +i 

ti(n) = Ski-. -Ps / Pt n {a l (t) - a l (t n )) flj(f) dt. 



III. 4 Quantum Langevin equations 

In this article what we call quantum Langevin equation is actually a restricted 
version of what is usually understood in the literature (cf [G-Z]); by this we mean 
that we here study the so-called quantum stochastic differential equations as de- 
fined by Hudson and Parthasarathy and heavily studied by further authors. 
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This type of quantum noise perturbation of the Schrddinger equation is ex- 
actly the type of equation which we will get as the continuous limit of our Hamil- 
tonian description of repeated quantum interactions. 

The aim of quantum stochastic differential equations is to study equations of 
the form 

dU t = J2 L)U t da){t), (10) 

i,j6AU{0} 

with initial condition Uq = I. The above equation has to be understood as an 
integral equation 

U t = I+ [ £ L)U t da)(t), 

J ° i,j€Au{0} 

for operators on Hq <E> <E>, the operators L % - being bounded operators on Ho alone 
which are ampliated to Ho <E> 

The main motivation and application of that kind of equation is that it gives 
an account of the interaction of the small system Ho with the bath $ in terms of 
quantum noise perturbation of a Schrodinger-like equation. Indeed, the first term 
of the equation 

dU t = L° U t dt + ... 

describes the induced dynamics on the small system, all the other terms are quan- 
tum noises terms. 

One of the main application of equations such as (10) is that they give explicit 
constructions of unitary dilations of semigroups of completely positive maps of 
B(Ho) (see [H-P]). Let us here only recall one of the main existence, uniqueness 
and boundedness theorems connected to equations of the form (10). The literature 
is huge about those equations; we refer to [Par] for the result we mention here. 

Theorem 12. -If Ho is separable and 

\\l\\ = ( yj ||l;h 2 ) 1/2 < +oo, 

;,j€Au{0} 

then the quantum stochastic differential equation 

U t = I + V] I L)U t da){t) 

admits a unique solution defined on the space of coherent vectors. 

The solution (Ut) t>0 is made of unitary operators if and only if there exist, 
on Ho, a self-adjoint operator H, operators Li, i e A and operators Sj, i,j G A 
such that (Sj)ij e A is unitary and the coefficients L l - are of the form 

L° = -(iH+^J2 L *M 

keA 
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keA 



L) = S)-5 ij I. 



IV Convergence theorems 

IV. 1 Convergence to quantum Langevin equations 

We are now ready to assemble together all the pieces of the puzzle and prove 
that the Hamiltonian dynamic associated to repeated quantum interactions spon- 
taneously converges to a quantum Langevin equation under some normalization 
conditions on the Hamiltonian. Notice that we no longer assume that Ej{h) has 
been conveniently constructed for our needs; in particular K, is not assumed to be 
unitary. 

Let /ibea parameter in M + , which is thought of as representing a small time 
interval. Let 2L(/i) be an operator on Ti <S>Ti., with coefficients lL l j(h) as a matrix 
of operators on Hq. Let u n (h) be the associated solution of 

u n +i(h) = E n+1 (h)u n (h) 

with the same notation as in section II. 3. In the following we will drop dependency 
in h and write simply iL, or u n . Besides, we denote 



1 



for all i,j in A U {0}. That is, 



£io — eoj — £ij — 0, £oo — 1- 

Note that from now on we take the embedding of T$ in $ for granted and we 
consider, without mentionning it, all the repeated quantum interactions to happen 
in Tl>(h), the subspace of $ associated to the partition S = {U = ih; i G IN}. We 
also make the convention that the default summation sets for sums is iLU {0}, e.g. 

Theorem 13.- Assume that there exist bounded operators L 1 -, i,j e A U {0} on 
Ho such that 

J2\\ l j\\ 2 < +o ° 



i,3 

and 



hm > 

i,3 



h £i i 3 



= 0. 



for alli,j = 0, . . . , N—l. Assume that the quantum stochastic differential equation 

dU t = J2 L Wtda)(t) 
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with initial condition Uq = I admits a unique solution (Ut) t>0 which is a process 
of bounded operators with locally uniform norm bound. 

Then, for almost all t, for every </>, tp in L°°([0, £]), the quantity 

< a <g> e(4>) , Psu [t / h ]Ps b <g> e{ip) > 

converges to 

<a<S>e((f)) , U t b®e{^) > 

when h goes to 0. 

Moreover, the convergence is uniform for a, b in any bounded ball ofH, uni- 
form for t in a bounded interval of M + . 

Remarks 

- This is where we particularize the index zero : the above hypotheses of conver- 
gence simply mean that, among the coefficients of iL, 

(_ZLq(/i) — I)/h converges, 

converges if either i or j is zero, 

lL l j(h) — Si j converges if neither i nor j is zero. 

We here meet the announced three time scales appearing in the Hamiltonian. We 
shall discuss the physical meaning of these normalizations in next section. 

- In the case where the operator iL is unitary and satisfies the convergence as- 
sumptions of the above theorem, then one can see that the limiting operators L*- 
are of the form given in the second part of Theorem 1 1 . 

- In that case, the solution (U t )teR+ en j°ys a particular algebraic property which 
we won't define here: it is a cocycle (see [H-P] or chapter 6 of [Pa3]). This property 
traduces the fact that the evolution of the system is, in the limit, memory-less. 

Consider the quantum stochastic differential equation (E) on Ho <S> 

dU s = J2 L Wsda)(s) 

where the U- are the bounded operators on the initial space Hq given by our 
assumptions. 

We consider that h is fixed and the associated partition S = {0 = t < ti = 
h < . . . < tk = kh < . . .} is also fixed. Note that we have chosen a regular 
partition only for simplicity and that all our results hold with general partitions 
when the mesh size tends to 0. We fix some bounded interval [0, T] of M+. 

We will proceed by successive simplifications. Consider the operator on the 
atom chain defined by 

w k = P s U tk P s . 
The following lemma will be used over and again. 

Lemma 14. -For any r < s, any vectors a £g> e{4>), b <E> with 

4>,%j) e L 2 (1R + ;H') nL°°(lR + ;H'), 
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we have 

\<a® e(4>) , (Us -U r )b® e(ip) >\<C \\a\\ \\b\\ (s - r) 

where C depends only on \\L\\, defined in Theorem 12, and on the L 2 and L c 
norms of and ip. 

Proof 



|< a <S> e((j>) , (U s - U r )b ® e(ip) >\ 



<||L|| / ||0(u)|||^(u)||||a(8)e(0)||||^ tt 6(8)e(V')||du 

J r 

from which the estimate follows, using the fact that ||0(n)|| < H^H^ and ||^(u)|| < 
HV'lloo' ^ ne are bounded and U is locally uniformly bounded. ■ 

The following lemma shows that (wk) k converges to (Ut) t>0 , in the same weak 
sense as in the theorem, as h goes to 0. 

Lemma 15. -For any t k < s, any vectors a <E> £(4>)> b <E> e(V0 with 

e L 2 (M + )nL 00 (lR + ), 

we have 

\<a(g) e((p) , (tufc - C/a)6 <g> e(V>) >| 

< C \\a\\ \\b\\ ((s - t k ) + - P s )e(<p)\\ + - Ps)e(il,)\\) • 
where C depends only on ||L|| and on t/ie L 2 and L°° norms of 4> and ip. 

Proof 

We have 

| < a <g) £(0) , (w k -U s )b <g) e(V') > | < 

< |< a <g> £(0) , (C/" tfc - C/s)6 <g> e(V>) >| 

+ \<a®e(<i)), (P s U tk P s -U tk P s )b®e(ij)>\ 
+ |<a®£(0), (C/ tfc i 3 5-C/ tfc )6®e(V')>l 



< 53 / l&( w )| |<a®e(» , L)U u b ® e(t/)) >\ du 



+ \\(I-P s )a®e(<t>)\\ \\U tk P s b®em\ 
+ j|C/ t ;a®e(0)|| ||(J-P 5 )6®e^)|| 
and we conclude as in the previous lemma. 
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We can now prove Theorem 13. 
Proof of Theorem 13 

Let ujj(h) be such that 

IL){h)-5 l3 I = h^{L)+u)(h)) 
for all i,j in A U {0}. In particular we have that, 

i,jeAu{o} 

converges to when h tends to 0. 

Consider the solution {u n ) ne]N , of 

with the notations of section II.3. Note that if A denotes the matrix iL — (Sijl)ij 
we then have 

U n +1 ~U n = A n+1 U n . 

Let F be the matrix [h eii V- + S^hl,®)^ where 



_ Jl ifz = jand(z,j)^(0,0), 
°^ \0 ifz^jor (z,j) = (0,0) 



and consider the solution (v n ) ne]N of the equation 

V n +1 -V n = F n+1 V n . 

Note that 

A n+1 = J2 A j a )( n + 1 ) 

and similarly for F n +\. 

Also note that a l An + 1) commutes with u n (resp. v n ), for they do not act 
on the same part of the space T$. Thus we get another useful way to write the 
above equations in terms of the basis a* (n) : 

u n +i -u n = y^^tt n a)(n + 1). 

and 

v n+1 - v n = (h £lj L) + hjhL ^ v n a)(n + 1). 

From the above lemma it is enough to prove the convergence to zero of u n — w n . 
We actually start with w n — v n . 
From the fact that 

U tk+1 -U tk =J2 f k+1 L)U s da){s) 

and thanks to the formulas for projections of Fock space integrals onto the toy 
Fock space in Theorem 11, one obtains the following expression for Wk+i — Wk (be 
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careful that the da^t) integrals give rise to a\{k) terms for all i, for / = £V a\): 



w k+1 -w k = J2 h ^ L ) (\ P s I ^Pt.Ut dt) a){k + 1)P S 
i,j¥(o,o) \ n Jt k J 

+ hL °o (l p s I h+1 ^ dt ) a i( k + 

+ ££ P *(r r +1 (LiPt k U t (ai(t)-4(t k ))) dt 



1 



tk 



[L]P th Ut{£{t) - aj(f*))) * + ^ 



As a consequence 



fc<n (i,i)/(0,0) V ' ,/tfc 



+ aj(fc + 1) ( ^Ps / ?7t eft^s - "/ 

k<n i 



■k + l 



+ EEI> 

fc<n i6A j'6A 



L*P tfe ^K(t)-a^ fc ))) eft 
+ i jT* fc+1 (L;P tfc ^(of (t) - a?(f fc ))) eft) P5 a}(* + 1). 



(11) 



We first wish to replace ^ P5 J^ k+1 Pt k Ut dt or P5 J^ k+1 Ut dt in the first two terms 
by wjfe. Lemma 14 allows us to estimate the error term. Consider two essentially 
bounded functions 0, •0 in £ 2 and two vectors a, 6 in the unit ball Hi of H; 

we expand the first terms in (11) and estimate each term : 



k<n 



<J2\Uk)\\Mk) 

k<n 



W-^-(* + i)^W> 

2 ptk + 1 



<a®e(0 fc ), Psh^'L)-} (U s - U tk )dsb ® e(tp k ) > 



where we have omitted a uniformly bounded factor 

a®e(^ fc+2 ),6®e(j/i[ H 2) 



From the fact that 



< Vh 



00' 



estimation of the error term of the form 

C llall 11611 Vh 



< Vh 



we obtain an 



for some constant C which depends only on the ||L|| and on the L 2 and L°° norms 
of 4> and ift. 
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On the other hand, the error term obtained by replacing ^ P$ J^ k+1 U t alt by 
Wk in (11) is clearly dominated by C h in norm because £V hL^aKk + 1) is just 

We now seek to evaluate the third sum in (11); for that consider again two 
functions 0, ip in L°°([0, t]) and two vectors a, b in the unit ball Hi of 7Y; we have, 
up to a uniformly bounded factor, 



£ <a® £(0) , iz7 P 5 ( / fc+ P tfc ^(ao(«) " a o(^)) 4^ + X ) Ps b ® e ^ 

i,jeA " Jtk 



> 



I ~ ~ ~ 1 /"*fc + l . „ 

= 2^ minify <a®e(<f> k ), -L l U s (a J (s) -a J (t k ))b®e(i(j k )ds> 

<Ch^\\L\\ 

for Oq(s) — a J (tk) is bounded on $ tfc , with norm y/s — tk- One obtains similarly 

>t k 



<a® e{4>) , T L° P s 
a 



■k + l 



PtMaP^s) - a^tk)) ds aUk + 1) P s b ® e(V>) > 



oo ' 



E 

ijeA 

<h" 2 \\L\\ Ul 

so that the third sum in (11) is bounded by C Vh2t ||0|| 

We have shown that, putting F k = E+ j{h £lj L 1 - + d^hL^) a l j(k + 1), 
|< a <g> e(0) , Ps(«; n - v n )P s b <g> e(V>) >| 



+ o(l) (12) 



= ^ < a ® e(0) , F fe (w fc -Vk)b® e(V>) > 

where o(l) is a term which converges to zero as h goes to zero uniformly for a, b 
in Hi and for t in a bounded interval. That uniform convergence property will be 
important in the sequel. 

Expanding Fk in equation (12) gives 



< a <g) e(4> n ) , (w n - v n )P s b <g> e(ip n ) > 

>j(k) ijji(k) < (L))* a ® e{<t) k ) , {w k - v k )b® e{ip k )> 



< E E 

i,j k<Cn 

+ J2J2 h \&W <(^o)*«®e(^), (w fc -v fc )6®e(^ fc )> 

where in each term we have omitted an uniformly bounded factor and the notation 
o(l) indicates a function of /j which converges to zero as h goes to zero. 

Since and ^ are essentially bounded, the quantities J+ ■ /i £iJ <f>j(k) ipi(k) 



and EieA^ &(*0 ^i( k ) 



are again of order h at most, with an estimate which is 
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independent on k. Besides, remark that normalizing all operators would only imply 
an additional constant factor, so that we can 
In that case the above implies that for all n, 



an additional constant factor, so that we can assume all the L*- to be contractions. 



SU Pa,beHi 



<a<g> e(4> n ) , (w n - v n )b <g> e{ip n ) > 

<a<S> e{4>k) , (w k - v k )b <g> e(tf) k ) > 



+ o(l) 



< frC^sup a6eWl 

k<n 

for some constant C. Here we have used our earlier remark that all convergences 
are uniform in a, b e Hi. This implies that 



su Pa,6€Wi 



< a <g> e(<j> n ) , (w n - v n )b <g> e(^ n ) > < (1 + Ch) n x o(l) 



and since nh converges to t, the quantity (1 + Ch) n is bounded so that 



sup a , b6Wl 



< a <g> e(<fi n ) , (w n - v n )b ® e(V> n ) > 



converges to zero as h goes to zero. 

We have proved the desired convergence property for the process (w k ) k >o. 
Now we will prove that 



sup„, b6Wl 



< a <g> e(4> n ) , (u n - v n )b ® e(^ n ) > 



converges to zero as h goes to zero. We have 



= F k{u k -Vk) + J2\Y1 hSxj H( h ) - h5,jL° )a)(k + 1) | a, 

k<n \ i,j 



k<n 



\F k + Y,h eii (uj(h)-h8ijLo)a.)(k+l)\ (u k -v k ) 



k<n 



^h^{u){h)-h8 i:i Ll)^{k+l) | v k 



(13) 



fc<n \ i,j 



Replace for simplicity ujj(h) by ujj(h) — HS^Vq (which does not change the validity 
of the earlier estimates) . The interest of the second form lies therein, that the term 
without recurring u k — v k can now be estimated thanks to our previous result: 

< a <g) e(0 fc ) , ^2 XI h £ii u)(h)a){k + 1) \v k b <g> e(^ fc ) > 

fc<n y i,j J 

is bounded by 

J2 w2 h£i3 \^( k )\ \^ k )\ \H( h )\\ su Pa,6e«i |<a®e(0 fc ), v fc 6®e(^ fc )> 

fc<n y i,j / 

(14) 
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with 

SUP a , 66Hl 



< a <g) e(0 fc ) , i> fc 6 <g> e(^fc) > 



< a <g> e(0fc) , (ufc - w k )b <g> e(ip k ) > 



+ 



< a <g> e(0 fc ) , w k b <S> e{ip k ) > ) 



and our estimate shows that the first term on the right-hand-side converges to 
zero uniformly in k as h goes to zero. The second term is, in turn, bounded by 
||£(V')|| since any w k is PsUt k Ps and as such has bounded norm. 

Thanks to our assumptions on the perturbative operators, the bound (14) we 
are interested in converges to zero as h goes to zero, uniformly for a, b in Hi- 

Besides, for the recurring term in (13) we obtain as before 



SU Po,6€Wi 



< a ' 



> 



(w k - v k )b<S> e(V> n ) 
< a <g> e(0fe) , (vk - w k )b ® e(^ fc ) > 



i,3 



k<n 

<hCj2 SU Va,b 



6H1 



thanks to the fact that the operators u e are assumed to have norms which converge 
to zero uniformly with h. We conclude as in the previous case. 

This ends the proof. ■ 

Under some additional assumptions, which are verified in many applications, 
we can very much improve the convergence. 

Theorem 16.- Consider the same assumptions and the same notations as in The- 
orem 13. If furthermore \\uk\\ is locally uniformly bounded, thenu^m converges 
weakly to Ut on all Ho <E> 

Proof 

Theorem 13 allows us to perform a e/3 argument with an approximation of 
any vectors of Ho®T3? by combinations of vectors a®e(4>), b<S>s(ip) with essentially 
bounded functions <j), ip. I 

One of the main application of this last theorem is the case when the matrices 
Ej(h) give rise at the limit to a matrix L such as in Theorem 12 (the case of a 
unitary solution (U t ) t>0 ). We shall show that the associated discrete evolution 
(u n ) ne ]N satisfy the conditions of the above theorem, so that the convergence of 
u [t/h] towards U t is weak. 

Theorem 17. - Consider a matrix E, on Ho <8> H with coefficients 

!L G = I-h{iH+\Y,LlL k ) + h^ 



= VhL, + hu° 
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K = -^hJ2 L l S J+ huJ o 

k 

Eij = I + Si- — SijI + huj 

where H is a bounded self-adjoint operator, (Sj)ij e A, is unitary, the Li, i e A are 
operators on TCq such that L* k Lk converges strongly, the coefficients ujj are such 

that Ylii j H^] || i s uniformly bounded and \\uj\\ Q {h) converges to when h tends 
to 0. 

Then the solution (u n ) ne]N of 

is made of invertible operators which are locally uniformly bounded in norm. 

In particular u^t/}^ converges weakly to the solution U t of the quantum stochas- 
tic differential equation (E). 

Proof 

A straightforward computation shows the special form of E induces many 
cancellations when computing the coefficients of E* E — I and of EjEj* — 7, and 
that they are of order h. Thus for h small enough the operators Ej* Ej, EjEj* and 
thus E are invertible. Thus so are the operators u n . 

Furthermore the above estimates show that ||iL|| < y/l + Ch. This easily 
gives the locally uniform boundedness of (u n ) ne]N and thus the desired weak con- 
vergence. ■ 

Specializing some more will allow us to answer the natural question of con- 
vergence of Heisenberg evolutions of observables: 

Corollary 18.-7/ the operator E is unitary and satisfies the conditions of Theo- 
rem 13 then the solution (Ut)teR+ °f (E) is unitary. In this case the convergence 
°f u [t/h] t° U t is strong and for all bounded operator X on Ho <E> almost all t, 
the sequence u* t ^ h ^Xu[ t /f^ converges weakly to U^XUt as h tends to zero. 

Proof 

It is easy to see from the above conditions that the operators (Lj)ij satisfy 
for all i, j 

L ] + L Y + L k L k* = 

fceA 

L) + L)* + J2 L * L ^ =° 

which implies that the equation (E) is of the form which has unitary solutions (see 
Theorem 12). 

By Theorem 17 and Proposition 1, for almost all £, u t / h is a sequence of 
unitary operators that converges weakly to a unitary operator, so that strong 
convergence also holds. 

It is now straightforward to prove the statement regarding convergence of 
u? t , h -,Xu[ t /h] to U^XUt for bounded X. ■ 
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IV. 2 Typical Hamiltonian: weak coupling and low density 

We are now coming back to the initial physical motivations of Theorems 13, 
16 and 17. These theorems show up quite strong conditions on the unitary op- 
erator Ej = e~ thH and a natural question now is: what kind of Hamiltonian H 
will produce such conditions on iL? What is the typical Hamiltonian for repeated 
quantum interactions which will produce quantum Langevin equations at the con- 
tinuous limit? 

In this section we answer partly that question. We answer it as we exhibit 
a large family of such Hamiltonians and we conjecture that they are the typical 
ones. We do not fully answer the question for we are not able to prove that they 
are the only ones. 

We keep here the notations of section II.1-II.3. 
On Hq®H consider the following Hamiltonian 

H = H <g> / + / ® H s + 4= ( Vi ® a * + V i ® a o) + \ Yl Di i ® a ) ( 15 ) 

^ h teA h i,jEA 

where Hq, Vi and are bounded operators on Ho (with Hq hermitian and 
Dij = D^), and Hs is bounded hermitian on H. 

The condition that H$ is bounded can be felt as a weakness in our conditions. 
But one has to keep in mind that 7i is just a "small piece" of the bath system; it 
is only ®jr+7Y which represents the bath. In general H is finite dimensional and 
the resulting continuous field is a Fock space. 

Put H! to be the closed subspace of TL generated by the basis elements X\ 
i E A; that is, the orthogonal of X° = O. Consider the "column operator" 




as an operator from Hq to Ho ® H' . Assume that this operator is bounded. 
The adjoint of V is then the "row operator" 

V* = (V 1 V 2 ...) 

from Ho®H' to H . 

Define the "matrix operator" D = (Dj) as an operator from Ho <8> H' to 
Ho <E> H' . We also assume D to be bounded. 

Let us relate the above Hamiltonian with the usual literature on weak coupling 
and low density limits. 

Recall that, in the literature about weak coupling limit, the bath is usually 
made of several harmonic oscillators (a Fock space) with associated creation oper- 
ators a*(g) and annihilation operators a(g), where g runs over a Hilbert space H. 
The Hamiltonian which is considered is then of the form 

H = H ® I + I ® H s + A (V <g) a*{g) + V* <8> a{g)) 
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The part V <S> o*{g) + V* <8> a(g) corresponds to the typical dipole Hamiltonian 
usually considered in the weak coupling limit or van Hove limit (cf [Dav], [D-J]). 

This Hamiltonian meets (15) when considering an ortho normal basis fa) of 
7i and when one writes 

(V <g> a*(g) + V* <g> a(g)) = J2(<^,g>V® a* fa) + < g , e z >V* <g> afa)) 

i 

Our time renormalization term 1 / \fh corresponds to the usual time renormalisa- 
tion for the weak coupling limit (t/X 2 = r). 

On the other hand, interaction Hamiltonians of the form D ® a*(f)a(h) + 
D* <g> a*(h)a(f) are typical of the low density limit (cf [APV], [AFL]). They meet 
(15) when decomposing in the same way as above: 

(D <g) a* (f)a(h) + D* <g> a* (h)a(f)) = ^(<e,, f><h, e$ >D ® a*fa)a(ej) + 

+ < f , e< >< ej , h >D* ® a*fa)afa)) . 
Our time renormalization 1/h is also the typical one for this limit. 

We are now back to our general Hamiltonian (15). 

We shall take compact notations for the unitary quantum Langevin equations 
of Theorem 12. Consider the equation 

dU t = J2 L Wtda)(t) (16) 

with 

L° = -(iK+\j2 L k L k) 



2 

keA 



T u — T 



^o = -E lis 1 : 

keA 

L l j = Sj — SijI, 

where K is a bounded self-adjoint operator and (S^)i,jeA is unitary. 
We write W for the column operator 



L 



2 



and S for the matrix operator [S]) . .. Then, with obvious notations 

L° = -(iK+ )-W*W) 



L° = W 
L = -W*S 
L=S-I. 
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With the Hamiltonian H given by (15), put 

U = e~ ihH 

and consider the evolution equation for repeated interactions associated to U : 

Un+l = U n+ lU n . 

Theorem 19.- The solution {u n ) n£]N * of the discrete time evolution equation 
converges strongly in $ to the solution (U t ) t>0 of the quantum Langevin equation 

dU t = J2 L Wtda)(t) 

where, with the same notations as above 

K = H + <Q, H s tt>I + V*D- 2 (sinD - D)V 
W = D- 1 (e~ iD -I)V 
S = e~ iD . 

Moreover, for any bounded operator X on Ho <8> 3>, u* n Xu n converges to U^XU t . 

The convergences are meant, as in Corollary 18, for almost all t. The ex- 
pressions D~ 2 (sinD — D) and D~ 1 (e lD — I) have to be understood has a short 
notation for the associated convergent series, even if D is not invertible. 

Proof 

Put kj = <X j , H$X l > for all i,j G A U {0}. We consider the column 
operator 

\k) = 



the row operator 

and the matrix operator 



k%I 
\ ) 

{k\ = {kll fcgj ...) 



k = (k)I). . . . 



H = 



They all are bounded operators. 

Put H = H + fcg/ and M = (Mj) . . with M\ = S^Hq + kU. We then have 

^+|*> ±D + M 

as an operator on Ho <8> W which is decomposed as an operator on 7i <8> (<Cf2© W). 
In particular 

hH VhV* + h(k 

VhV + h\k) D + hM 
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Let a be a bound for 



H 



\V\\, \\D\\, ||M||, ||<fc|||, |p)||. 



Lemma 20. -For all m e IN we have 

hA m + h 3 / 2 R m VhB m + hR 2 m 

^hC m + hRl D m + hRi 



(hH)' 
with 



m 



for all X = A,B,C,D and 



\XrnW < a 



\Kn\\ < 7 m_1 « m 



for all i = 1, 2, 3, 4; with 

A m+1 = V*C m7 A = I 7 A 1= H 

B m+1 = V*D m , B = 0,B 1 = V* 

Cm+i = DC m , Cq = 0, C\ = V 

D m+1 = DD m , D = I. 

Proof of the lemma 

For m = the statements are clearly satisfied. For m = 1 we find the 
announced A 1 ,B 1 ,Ci,D 1 and R\ = 0, R\ = (k \, Rf = | k), Rj = M. The norm 
estimates are then clearly satisfied. Now, applying the induction hypothesis and 
computing (hH) rn+1 from (hH) m we get 

A m+l = V*C m and \\A m+1 1| < aa m = a m+1 

R m+1 = VhHA m + hHR m + V*R m + (k \C m + Vh(k \R m 

\\R m +i\\ < 5(7 m " 1 a m+1 ) < 7 m a m+1 

B m+1 = V*D m and ||S m+ i|| < aa m = a m+1 

R m+1 = VhHB m + hHR 2 m + VhV*R m + (k \D m + h(k \R m 

\\R m+ i\\ < 5(7 m " 1 a m+1 ) < 7 m a m+1 

C m+1 = DC m and ||C m+ i|| < aa m = a m+1 

R m+1 = \fhVA m + hVR m + DR m + VhMC m + hMR m 

\\R m +i\\ < 5(7 m " 1 a m+1 ) < 7 m a m+1 
D m+1 = DD m and \\D m+1 \\ < aa m = a m+1 

R m+1 = VB m + hVR 2 m + Vh\ k)B m + h\ k)R 2 m + DR m + MD m + hMR m 
\\R A m+i\\ < 7(7 m " 1 a m+1 ) < 7 m a m+1 . 
The lemma is proved. 
In particular, we get 





= D m 


B m 


= V*D rn ~ 1 


C m 


= D m ~ 1 V 


A m 


= V*D m ~ 2 



V. 
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By the norm estimates of Lemma 20 we have that the series 

(-,)•» 



E 



-RL 



ml 



are norm convergent. Let us denote by R k their repective limit. 
We get 

(-*) 



ml 



■{hH)" L = 



V 

This gives 



+h s / 2 R 1 

\oo ( — i) r 
IT 

+hR 3 



+hR 2 



oo (-i) r 



I + v°° 



/ 



/ / — ihH + hV*D~ 2 (e~ iD — I + iD)V VhV*D- l (e~ lD - I) \ 
+h 3 / 2 R 1 +hR 2 



U 



JhD- x {er iD - I)V 



\ 



■l(„-iD 

+hR 3 



I+(e-* D -I) 



+hRi 



I 



We are exactly in the conditions of Corollary 18 and we thus get the strong con- 
vergence to the solution of 

alU t = J2 L Wtda)(t) 



where 



7"° 
L° 

L 



-iH + V*D- 2 (e~ iD -1 + iD)V 
D-\e- lD -I)V 
V*D- 1 {e~ iD - I) 
e~ iD - I. 

If we put W = L° and S = e~ lD , we then get 

-W*S = -V*(e lD - ^D^S = V*D-\e- iD — I) = L 

and 

--W*W = --V*(e lD - I)D- 1 D- 1 (e~' D - I)V 
2 2 



= -)-V*D- 2 {cosD -I)V. 



This shows that 



1 



L° Q = -iH - iV*D~ 2 (smD - D)V - -W*W. 



The theorem is proved. 
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Let us interpret the above theorem in terms of weak coupling and low density 
limit again. If in the above Hamiltonian we consider no term Dij, that is, 

H = H <g) I + I <g> H s + T] (Vi <g> a° + V* <g> aj) 

then we are in the usual situation of a weak coupling limit, with its typical dipole 
interaction Hamiltonian. The quantum Langevin equation we obtain in Theorem 
16 then simplifies to 

dU t = - (^H + K O , H s n >I+lYl V i V ^j Ut dt+ 

+ J2 V*U t cfaj(f) - J2 V * U t da (t). 

i i 

This is also typical of the "diffusion" type of quantum Langevin equation one 
can meet in the literature for this kind of limit. The fact that only creation 
and annihilation quantum noises are involved is here the quantum analogue of a 
classical stochastic differential equation with Brownian noise. Note that if the Vj's 
are such that V* = —Vi then the above quantum Langevin equation becomes a 
classical stochastic differential equation with Brownian noises: 



dU t = 



We refer to [Atl] for a complete discussion on the classical stochastic interpreta- 
tions of the quantum noises. 

On the other hand, if in the Hamiltonian we consider no term Vi, that is, 



1 \ ^ 

H = H ®I + I®H S + - D i: j® a 



h ^ ~ " J ~ ^ J 

i,j€A 

then we are in the usual situation of a low density limit, with its typical scattering- 
type interaction Hamiltonian. The quantum Langevin equation we obtain in The- 
orem 18 then simplifies to 

dU t = -i(H + < O , H S Q >I)U t dt+ ( S j ~ 5 v) U t da){t) 

ijeA 

where S = e~ tD . 

This is also typical of the "Poisson" type of quantum Langevin equation one 
can meet in the literature for this kind of limit. The fact that only "exchange" 
quantum noises are involved is here the quantum analogue of a classical stochastic 
differential equation with Poisson noises. 

Our formalism enables us to handle these two types of limits in a single 
setup, and this could not have been done at the classical stochastic calculus level. 
This seems to be the first time in the literature, for, to our knowledge, weak 
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coupling limits and low density limits have always been considered separately, as 
very different objects. 

The surprise in our setup in the apparition of the term 

V*D~ 2 (smD - D)V 

only when both the limits are in presence in the Hamiltonian. 
Indeed, the term 

L° Q = -(iK+ ^W*W) 

is the driving term of the dynamic associated to the quantum Langevin equation. 
In some sense it is the generator of the dynamic on Ho- The part ^W*W is 
representative of the dissipation from Ho to the bath. The part iK is an effective 
Hamiltonian on Hq. The apparition of this new contribution V*D~ 2 (sinD — D)V 
is new, and we have no physical interpretation of it. 

It is just clear that it results from the combined effects of the weak coupling 
limit and the low density limit. 

IV. 3 Hamiltonian description of von Neumann measurements 

Our setup and approach allows to construct an Hamiltonian description for 
the usual von Neumann measurement procedure (collapse of the wave packet pos- 
tulate). 

On some quantum system state space Ho consider an observable A with dis- 
crete spectrum (maybe infinite). Let Pi, P 2 > ••• denote the associated spectral (or- 
thogonal) projections, with J2k P k = I. 

We want to give a model for the action on Ho of an exterior measurement ap- 
paratus which measures the observable A. That is, the action of the measurement 
apparatus on Ho should be to transform any observable X of Ho into 

PkXPk. 

Let H be a Hilbert space with one more dimension than the number of projec- 
tors Pk involved above (infinite dimensional if the IV s are in infinite number). On 
H consider an orthonormal basis O = X , Xi, X 2 , ... and the associated creation 
operators a° and annihilation operators a§, as in sections II.1-II.3. Consider the 
following Hamiltonian on Hq <S> H: 




Let U = e lhH and consider the process 

of repeated quantum interactions on Ho <8> 0iv« 

Theorem 21.- The repeated quantum interaction process (u n ) ne]N * converges 
strongly inHo®<S) ]R + when h tends to 0, to the solution (U t ) t>0 of the quantum 
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Langevin equation 

alU t = -^U t dt + J2 {i p kU t da° k (t) - iP k U t dag(t)) . 

k 

Furthermore, for any observable X on Tio, the partial trace along the field T~l, 
in the vacuum state, ofU^(X <g> I)Ut converges, when t tends to +oo ; to 

PkXPh- 



Proof 

In the basis O = Xq, Xl, X 2 , ... for H the Hamiltonian H writes as 



H 



In particular, as an easy direct computation shows, we have 
/ cos y/h I — sin y/h P\ — sin y/h P2 
sin y/hPi cos y/hP x 
sin y/h P 2 cos y/h P 2 



( 




~% p > 


•A 






















\ i 






'■•/ 



TL = 



V 



/ 



This unitary matrix clearly satisfies the conditions of Corollary 18 (one could also 
have directly applied Theorem 19 to the Hamiltonian), we thus have the strong 
convergence to the solution of the announced quantum Langevin equation 

dU t = -\u t dt + Y J { lP k U t da° k (t) - iP k U t da k (t)) , 
1 k 

with Uq = I. 

Now, consider the evolution under (U t ) t>0 of a system observable X: 

U* t {X®I)U t 

and the partial trace along the field <S)jr+ ^ m the vacuum state: 

p t (x) = <n,u t *(x® i)u t n > 

with the notation 

< a , < O , AQ > b > Ho = <a®0,A(6®0) > Wo ®g) n 

for any operator A on Ho <8> <S>j?+ an y a i b e ^o- 

It is well known form the usual theory of quantum stochastic differential 
equations (cf [H-P]), that (Pt) t >o is then a semigroup, on B(Hq), of completely 
positive maps with Lindblad generator 
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One could also have obtained this "quantum master equation" by simply com- 
puting the discrete time quantum master equation associated to («Ti) n6iV ., using 
Theorem 2 and then by passing to the limit h — > to recover (P t ) t>0 (see section 
IV. 5). Our approach has the advantage to also describe the exact equation for the 
interaction with the bath. 
We thus get 

k 

so that 

C 2 {X) = -C(X) 

and 

P t (X) = e tc (X) = (I + (1 - e-*)£) (X). 

That is, 

P t (X) = (1 - e"*) J2 p kXPk + e~*X 
k 

which converges to 

k 

when t tends to +oo. ■ 

We thus have proved the a von Neumann measurement apparatus can be 
described in an Hamiltonian setup by, first considering an Hamiltonian description 
of a repeated quantum interaction, secondly passing to the limit to continuous 
quantum interactions (h — > 0) and thirdly passing to the limit to large times 
(t -> +oo). 

IV. 4 One example 

We shall here follow a very basic example. It is actually the simplest non- 
trivial physical example and it already gives very interesting consequences. 

Assume Ho = 7~t = C 2 that is, both are two-level systems with basis states O 
(the fundamental state) and X (the excited state). 

During the small amount of time h the two systems are in contact and they 
evolve in the following way: 

- if the states of the two systems are the same (both fundamental or both 
excited) then nothing happens; 

- if they are different (one fundamental and the other one excited) then they 
can either be exchanged or stay as they are. 

In the basis {O <g> O, O <g> X, X <g> O, X <g> X} the operator iL is taken to be of 
the form 



iL = 



Z 1 











COS (X 


— sin a 





sin a 


cos a 


Vo 
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The associated Hamiltonian is thus 

/0 

H = 






Vo 





ia/h 






-ia/h 








0/ 



,-ihH 



so that E — e 

For the choice a = Vh, that is, 



H 



1 

7^ 






Vo 






i 




we get 



TL = 






Vo 





cos 

sin Vh 






— z 










0/ 





— sin v^/i 
cos Vh 
1/ 



Repeating this interaction leads to an Hamiltonian equation of the form 

U n +1 = lL n+ iU n 



with coefficients 





cos Vh 




I - sin Vh 



JL\ = 
, E\ 



sin Vh 


cos Vh 




for iL. We then have 



iL° 
lim — =L = 

h-o Vh 

lim Ej\-I = 

h— 





-1/2 

1 










We are exactly in the condition for applying Theorem 15 since for all h the matrix 
Ej{h) is unitary. We get that for almost all t, U[ t / h ] converges strongly to U t where 
(Ut)teR+ is the unitary solution of 



dU t = --V*VU t dt + VU t da^t) - V*U t dal(t) 



with V = 



. This equation is the well-known quantum Langevin equation 



associated to the spontaneous decay into the ground state in the Wigner-Weisskopf 
model for the two-level atom (see [M-R] for example). 
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IV.5 From completely positive maps to Lindbladians 

Recall, from section II. 3 that the solution (u n ) ne]N of the equation 

Un+l = JLn+lUn 

with uq = I induces a completely positive evolution on the small system. Namely, 
in the Heisenberg picture, one has for any a, b in Ho, any X in B(Ho), 

< a <g) ft , u* n Xu n b <g> O > = < a , £ n (X) b > 

where 

Theorem 22. -Let Ki(k) = [lL^{K)) i Q n be a family of matrices such that 
J2ilLi(h)*lLi(h) = I and such that 

- (K, (h) — I)/h converges to some Lq, 

- E^/Vh converges to some L° for all i = 1, . . . , n. 

Then there exists a self-adjoint operator H in 13 (Ho) such that for all t G M + , 

£[t/h] > e tc 

in operator norm, where C is the Lindblad generator 

C(X) = l [H,X]+ 1 -Y: (urXL" - LfL°X - XLfLi) 

and I is defined above. 
Proof 

It is a straightforward computation that 

n 

£(X) =X + h(L° *X + XL° ) + Lfxtf + o(h \\X\\). 

i=i 

The equality £(I) = I entails 
so that 

V ieiN J 

is self-adjoint. We denote it by H. Then £ is of the form 

£(X) =X + h (i[H, X] + ±J2 fa? XL* - L?L*X - XL? L°)\ + o(h \\X\\), 

which is, with the notations of the statements 

£(X) = X + hC(X) + o(h \\X\\). 
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The above mentioned convergence is therefore clear. ■ 

As a consequence, it is very easy to obtain approximations of solutions of 
continuous-time master equations: 

by solutions of discrete-time ones: 

Yet notice that the master equation gives no information whatsoever on the in- 
teraction between the small system and the environment or on the environment 
itself. On the other hand, the associated quantum Langevin equation contains the 
information of the whole system; this justifies our effort. 

Notice that, in the above proposition, no hypothesis is needed on the other 
coefficients of the matrix K,(h). Their properties actually depend on the choice of 
additional features of the matrices IL(h), for example their unitarity. The possibil- 
ity of choosing K/(h) to be unitary and obtain in the end the desired Lindbladian 
L is described by Parthasarathy in exercises 29.12 and 29.13 of [Par]. 

What's more, these manipulations show that the hypotheses of convergence 
of Theorem 13 are not as artificial as it seems, and are not only convenient as- 
sumptions we set up in order to obtain the right convergence. Indeed, to TL{K) is 
associated both a dynamic on the observables, as we have seen, and an evolution 
r, defined by 



< a , r n b > = <aCg>fi, , u n 6(g>fi> 
for all a, b in Ho, which turns our to be 

If one assumes that T^/ft] converges for almost all t and that -ZLq is assumed to be 
continuous at h = then the assumption on iLp in Theorem 16 is to be fulfilled; this 
implies that J2ieA-^o* -^o = ~H^o* + ^o) + °(^)> so that the other assumptions 
of convergence of Theorem 13 are natural. 

The other conditions described in Theorem 17 are in turn necessary if one 
wants the process (Ut) t>0 obtained in the limit to be unitary or alternatively the 
matrices _ZL(/i) to be sufficiently close to unitarity. 
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